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Deep-seafloor organisms consume oxygen, which can be measured by 

in situ benthic chamber experiments. Here we report such experiments 

at the polymetallic nodule-covered abyssal seafloor in the Pacific Ocean 

in which oxygen increased over two days to more than three times the 
background concentration, which from ex situ incubations we attribute to 
the polymetallic nodules. Given high voltage potentials (up to 0.95 V) on 
nodule surfaces, we hypothesize that seawater electrolysis may contribute 
to this dark oxygen production. 


Oxygen (O,) is prevalent in deep-sea surface sediments where its rate 
of consumption reflects the sum of aerobic respiration and oxidation 
of reduced inorganic compounds produced by anaerobic decay. These 
processes define sediment community O, consumption (SCOC), and 
quantifying SCOC is needed to estimate fluxes of major elemental cycles 
through marine systems! *. We undertook multiple in situ benthic cham- 
ber lander experiments to measure abyssal SCOC in the Nauru Ocean 
Resources Inc. (NORI)-D licence area of the Clarion—-Clipperton Zone 
(CCZ; Extended Data Fig. 1 and Extended Data Table 1) where polymetal- 
lic nodules cover extensive areas of seafloor. Sediments and nodules 
were exposed to different experimental treatments, which included the 
addition of dead-algal biomass, dissolved inorganic carbon and ammo- 
nium (NH,’) or cold filtered surface seawater. No-injection controls were 
also performed. Incontrast to previous deep-sea O, flux studies that only 
showed SCOC, we consistently found that more O, was accumulating in 
the chambers than was being consumed, resulting in net O, production. 

Constant linear decreases in O, optode readings were observed in 
two experiments (Fig. 1), and SCOC determined by in situ O, micropro- 
filing was 0.7 mmol O, m” dtindicating that SCOC occurs in NORI-D as 
in many abyssal habitats” *. However, O, concentrations in 25 benthic 
chamber incubations started at 185.2 + 2.9 pmol 1” (1 standard error 
(SE)) and reached O, maxima between 201 and 819 umol I over 47 h 
(Fig. 1), indicating net dark O, production (DOP) corresponding to 


rates of 1.7-18 mmol O, m” d”. Independent measurements of O, 
concentration using the Winkler method also showed DOP (Extended 
Data Fig. 2), providing evidence that the optodes were not malfunction- 
ing. Nostatistically significant difference in the total net O, produced 
(maximum [0,] - initial [O,]; Extended Data Table 2) was found between 
chambers (ANOVA, F, , = 0.107, p = 0.900) or experimental treatments 
(ANOVA, F; = 0.876, p = 0.489), ruling out any experimental bias. We 
found no difference in the total net O, produced between cruises 
(ANOVA, F, = 0.391, p = 0.684), though DOP was correlated to the 
average surface area of the nodules (Spearman’s correlation, p = 0.664, 
p=0.031). A re-evaluation of in situ O, optode data collected from 36-h 
benthic chamber experiments in the abyssal eastern and western CCZ 
(Extended Data Figs. 1 and 3) also showed DOP, indicating its occur- 
rence in multiple locations across the CCZ. Our findings contrast with 
all published deep-sea benthic O, flux studies and suggest that DOP may 
provide O, for benthic respiration. Whereas the DOP measured was 
greater than SCOC, we would urge caution when temporally upscaling 
our results, as the nonlinear production of O, suggests that DOP may 
not be continuous in nature. Moreover, the variance in DOP activity 
seen between experiments and its relationship to nodule surface area 
suggests DOP activity may change with nodule spatial density and type 
(for example, diagenetic versus hydrogenetic), so upscaling our results 
by area is also imprudent without additional studies. 
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Fig. 1| Oxygen concentrationsin pmol I measured by calibrated O, optodes 
through time in hin the different benthic chamber incubations. a-c, The in situ 
benthic chamber lander deployments were made during the 5D (a), 5E (b) and 7A 
(c) cruises to the NORI-D license area (Extended Data Fig. 1). Nodules were present 
inall incubation experiments. The green hue, blue hue and red lines in the 5D figure 
(a) denote dead-algal biomass, dissolved inorganic carbon + NH,’ and filtered 
seawater treatments, respectively. The gap in the optode data in AKS279-Ch.3 

was caused by the optode periodically not logging data. The black line indicates 
ambient O, concentration measured on the outside of the benthic chambers 
during AKS273 onthe 5D cruise. The green and yellow hue lines in the 5E (b) and 7A 
(c) figures denote the dead-algal biomass and control (no injection) treatments, 
respectively. The minor drops seenin some of the O, concentration profiles at 28, 
38 and 47 hare caused by the dilution of the chamber water with 50 ml of seawater 
that was entrained from the outside into the chamber through a 1.5 m (0.25 cm 
diameter) open tube when the syringe sampler collected seawater samples from 
within the chamber. The constant O, concentration measured during the first 2h 
of the 5D and 7A experiments was due to the stirrers being turned off for 1 h to allow 
the substrates (for example, dead-algal biomass) to sink to the sediment surface. 
Stirrers were turned on during the 5E expedition from the moment the lander was 
deployed until the lander returned and power to the stirrers was disconnected. 


Several lines of evidence indicate that the DOP was not caused by 
experimental artefacts. First, the total O, change between the experi- 
mental and control (non-injection) treatments was statistically indistin- 
guishable, and a steady increase in O, concentration was recorded over 
many hours in multiple experiments; these observations demonstrate 
that DOP was not attributable to the injection of exogenous fluids. 
Second, diffusion of O, from trapped air bubbles within the chamber 
was unlikely because each chamber uses two one-way valves in the lid 


to purge air fromthe chambers as the lander sinks. Even if an air bubble 
could be trapped long enough to reach the seafloor, gaseous diffusion 
of O, into the water phase would take <1 s at 4,000 m depth (Extended 
Data Table 3), whichis inconsistent with the steady increase in O, over 
many hours seen in multiple experiments (Fig. 1). Third, intrusion of O, 
from the plastic chambers into the water phase is unlikely (Methods) as 
they are built from polyoxymethylene, which is both highly inert and 
chemically stable in well-oxygenated settings and would not explain 
the variation in DOP because all experiments used identical materials. 
Last, DOP was also observed during 48-h ex situ sediment incubations 
(Extended Data Fig. 4). 

Several lines of enquiry were pursued to explain the DOP. Subsur- 
face advection of oxic bottom water from seamount flanks into seafloor 
sediments” and then into the chambers was discounted based on 
in situ O, microprofiling that showed pore water was a net sink for O, 
and undersaturated compared with the O, seen in the chambers. Fur- 
thermore, DOP was measured in sealed ex situ experiments (Extended 
Data Fig. 4) that prevented O, intrusion from below. It is unlikely that 
biological mechanisms were responsible for the bulk of the DOP as 
ex situ core incubations revealed DOP in the presence of poison (HgCl,; 
Extended Data Fig. 4). Whereas many microbes in the CCZ are able to 
detoxify Hg (II) to Hg (0)’, and some microhabitat pore spaces in the 
core may have remained HgCl, free, the taxa known to be capable of 
DOP (for example, Nitrosopumilus maritimus) are killed by its addition’. 
We also observed weak statistical support between the relative abun- 
dance of certain nitrogen-cycling microbial taxa and DOP (for example, 
Candidatus Nitrosopumilus p = 0.474, p = 0.420). The fact that DOP 
was detected in ex situ controls containing only polymetallic nodules 
(Extended Data Fig. 4) suggested that the DOP was linked to their pres- 
ence. Hence, we estimated the potential contribution of radiolytic O, 
production using a kinetic model’ and found 0.18 umol IO, would be 
generated by this process within 48 h. We also modelled the chemical 
reduction of manganese (IV) oxide at in situ temperature (1.6 °C) across 
a range of pH and O, conditions encountered at the seafloor to assess 
if this reaction (2MnO, > 2MnO + O,;Extended Data Fig. 5) could liber- 
ate the O, but found that <0.1 nmol of manganese (IV) oxide would be 
chemically reduced to manganese (II) at seafloor conditions. As such, 
localized radiolytic O, production from the sediments and nodules and 
chemical dissolution explain only a negligible proportion (< 0.5%) of 
the DOP observed. 

The oxygen evolution reaction requires an input voltage of 1.23 V 
plus an overpotential of approximately 0.37 V to split seawater into 
H, and O, (ref. 10) at NORI-D’s seafloor mean pH (7.41). This value can 
be lowered by several hundred millivolts if the reaction proceeds via 
the lattice-oxygen-mediated mechanism". Use of metal catalysts such 
as Mn oxides enriched with transition metals (for example, Ni) found 
in nodules” and characterized by large tunnel areas and abundant 
defect sites can optimize the adsorption of reactants and enhance 
conductivity and catalytic performance™®*, We tested the electrical 
potential between two platinum electrodes at 153 sites on the sur- 
faces of 12 nodules (Fig. 2) from the UK1, NORI-D and Bundesanstalt 
fiir Geowissenschaften und Rohstoffe (BGR) license areas. Although 
the potentials between different positions on the nodules were 
highly variable, - potentials up to 0.95 V were found and high mean 
background-corrected potentials were detected under cold-water 
conditions (Fig. 2 and Extended Data Table 4). On the basis of these 
studies and DOP being observed in nodule-only ex situ incubations 
(Extended Data Fig. 4), we hypothesize that the DOP may have partly 
resulted from seawater electrolysis, with the necessary energy coming 
from the potential difference between metal ions within the nodule 
layers, leading to an internal redistribution of electrons. Whereas 
questions remain concerning this potential mechanism (such as the 
identity of the energy source(s), longevity of DOP, catalytic stabilities, 
electrochemical conditions on exposed versus buried nodules surfaces 
and the influence of different chemistries within the nodule layers), 
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Fig. 2| Box and whisker plots of background-corrected voltage potentials on 
nodule surfaces. The nodules were collected from the NORI-D (1-5), UK1 (6-8) 
and the BGR (9-12) license areas. Potentials were measured at 21 °C (nodules 1-12) 
and 5 °C (nodules 6 and 7 cold) and between two different UK1 nodules (Tests 
land 2) and across the surface of a metamorphosed carbonate rock (control). 
Means are designated by the ‘x’ symbol, medians by the line, boxes show the lower 
and upper quartile values (excluding the median), whereas the whisker bars refer 
tothe minimum and maximum data values. The number of technical replicate 
measurements made at different points on the surface of each nodule/rock to 
make each box-whisker is shown by the number above each whisker bar. 
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the ‘geo-battery’ hypothesis was supported by the link between DOP 
and nodule average surface area. This connection could be due to an 
increased abundance of anode and cathode sites or a greater abun- 
dance of high Ni and Cu dendritic porous layers in larger nodules”. 
Assuming the ‘geo-battery’ is partly responsible for the DOP observed, 
the initial high DOP rate may have been related to the ‘bow-wave’ of the 
lander removing sediments from the surface of the nodules and expos- 
ing electrochemically active sites on the nodules. The slowdown in DOP 
seen later in the incubations could have then been caused by a reduc- 
tion in voltage potential and/or degradation of metal-oxide catalysts 
that has been observed in Mn oxide catalysts previously’. Whereas 
this process requires further investigation, if true, DOP activity may 
fluctuate with sediment coverage on the nodules inviting the urgent 
question of how sediment remobilization and distribution over large 
areas during deep-sea mining may influence DOP. 

Understanding the mechanism(s) behind DOP, its temporal nature 
and its spatial distribution will allow its role in abyssal ocean ecosys- 
tems to be better understood. Future studies of DOP in the deep sea 
may also shed light on broader relationships between metal-oxide 
deposition, biological evolution and the oxygenation of Earth’*”. 
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Methods 

A benthic chamber lander was deployed in the NORI-D license area six 
times in May-June 2021 (5D cruise), five times in November-December 
2021 (SE cruise) and five times in August-September 2022 (7A cruise) 
(Extended Data Fig. 1 and Extended Data Table 1). The lander comprised 
three independent, autonomous, square benthic chambers (484 cm’) 
separated by approximately <0.5 m. After arriving at the seafloor, 
the lander waited for 0.07-1.34 d before the chambers were pushed 
into the sediment to create an enclosed microcosm of the seafloor. 
Ten minutes into the incubation period, the enclosed chambers were 
injected with 50 ml of one of three solutions: (1) 0.45-m-filtered, cold 
surface seawater containing 79.2 mg of freeze-dried Phaeodactylum tri- 
cornutum algae, (2) 32 uM Na,HCO, and 40 uM NH, Cl dissolved in cold 
artificial seawater (salinity 35) and (3) 0.45-um-filtered, cold surface 
seawater. On some occasions, the injection mechanism failed allowing 
the response to control (no injection) conditions to be measured. The 
seafloor in the study area had a temperature of 1.6 °C + 0.006 °C (SE, 
n=28) anda pH of 7.41 + 0.05 (SE, n = 17). Immediately after the injec- 
tion, the overlying water was mixed with a submersible stirrer at 60 rpm 
for 1 min before the stirrer was turned off that allowed any particulate 
substrates to settle for 1h. After 1 h, the stirrer was then turned on again 
for the remainder of the experiments. During the SE expedition, the 
stirrers were programmed to continually stir the overlying water even 
immediately after injection. 

The syringe samplers removed approximately 50 ml of seawater 
from the water phase of each chamber at 0.1 or 0.03, 1, 3, 9, 28, 38 and 
47 h into the incubation experiment. Oxygen optodes (CONTROS 
HydroFlash O, manufactured by Kongsberg Maritime Contros GmbH) 
mounted in the lid of each chamber logged O, concentrations in the 
chamber every 10 seconds throughout each experiment. Two days 
before the first lander deployment of each cruise, the optodes under- 
went a two-point, multi-temperature calibration using 0 and 100% O, 
calibration solutions at 1.2, 7, 18 and 30 °C following the recommen- 
dations of Bittig et al. (ref. 18). On the 5D cruise, we also calibrated 
the sensors 2 d after the last lander experiment so we could estimate 
optode drift, which was negligible (0.27 umol Id“) over the course 
of the six-week cruise. The 0% and 100% O, saturation solutions were 
created by bubbling 0.45-m-filtered surface seawater ina bottle sitting 
ina water-chilling/heating unit with N, gas (0%) or an aquarium air bub- 
bling unit (100%) for 30 min. The O, concentration of the calibration 
solutions was confirmed in triplicate by Winkler titration. After incubat- 
ing seafloor sediments for 47 h, the lander chambers were closed by a 
shutter door at the base of the chambers, and the chambers were then 
pulled slowly out of the sediment, which took 1h. The lander was then 
recalled from the seafloor. In eight instances, the lander programme 
did not finish and the doors did not shut, preventing the sampling 
of sediment and determination of the volume of the water phase in 
the chambers (Extended Data Table 2). Once the lander was back and 
secured on deck, the chambers were opened and the water above the 
sediment removed via syphoning into a bucket. The distance fromthe 
top of the sediment to the base of the chamber lid was then measured in 
four places to get an accurate water depth for water volume estimates. 
Whenever possible, a photograph was then taken of the chamber sedi- 
ment and nodules from directly above the opening of the chamber. All 
syringes containing water samples were removed and taken to the 
shipboard lab for immediate processing or stored in a cold lab (4 °C) 
before processing. The optodes were removed and their onboard data 
downloaded to a computer. Finally, the nodules were removed from 
the chambers and washed of attached organic debris with cold (4 °C), 
0.45-um-filtered surface seawater and placed in sterile Whirlpak bags to 
be weighed in the laboratory later. The number of polymetallic nodules 
at the seafloor determined from chamber counts was 1170 + 97 m”. 

Unfiltered syringe sample seawater was carefully transferred 
from each 50 ml syringe to a 12 ml exetainer via a 10 cm tube attached 
to the syringe nozzle, ensuring no air bubbles were introduced and 


immediately fixed for microWinkler titration. The sample was then 
mixed thoroughly using a glass bead placed in the exetainer and placed 
inthe darkina4 °C refrigerator for 30-45 min to allow the precipitate to 
settle. Once the precipitate had sedimented, the exetainers were shaken 
again and left for 2-3 h before Winkler titrations were performed. All 
titrations were completed within 12 h after sampling to determine dis- 
solved O, concentrations. Each Winkler sample (approximately 5 ml) 
was titrated twice, and duplicate measurements showed minor differ- 
ences in O, concentration (5D cruise error: 3.5 + 0.3 pmol 17, n = 71; 5E 
cruise error: 1.3 + 0.2 pmol I", n = 69; 7A cruise error: 2.8 + 0.4 pmol I7, 
n= 84). Winkler O, concentration data were averaged for each syringe 
sample. The O, concentrations estimated by Winkler analysis were 
22 +1% (n= 42, SE, 5D cruise), 8+ 4% (n = 39, SE, 5E cruise) and 24 + 2% 
(n= 40, SE, 7A cruise) lower than the concentrations measured by the 
optodes at the same time point in the same incubations most likely 
due to out gassing of supersaturated O, caused by depressurization 
and warming of the externally mounted syringes (whose samples were 
used for Winkler analyses) during the lander recovery to the surface. 
Back on shore, the final O, concentration values were calculated 
following Bittig et al. (ref. 18) from the optode, calibration and in situ 
pressure data that was derived from the depth where each lander 
deployment was made. Time stamps in the optode data were compared 
to the lander computer programme times so the optode readings 
could be aligned to the schedule of the chamber experiment. The total 
change in O, concentration in each chamber was then calculated from 
the volume of the water phase above the sediment and the difference in 
O, concentration from when the chambers started to seal off the sedi- 
ment to the point when the maximum O, concentration was reached. 


Benthic O, microprofiling 

Benthic O, microprofiles were made during lander deployments 
AKS313, AKS316, AKS318 and AKS321 during the SE cruise using a 
UNISENSE deep-sea microprofiling unit mounted <0.5 m from the 
benthic chambers. The microprofiles were made using 20 cm O, 
microsensors that penetrated the sediment in 0.05 mm steps. The 
microsensors were calibrated 2 h before the lander deployments at 
in situ temperature (1.6 °C) at 0% and 100% O, saturation (above). At 
each sampling depth, the microsensor stopped for 5 s before each 
measurement was made. The sensor then recorded five individual O, 
concentration measurements. The average of these five measurements 
was taken for each depth point. The sediment surface was determined 
manually based on the turning point in the slope of O, concentration 
with depth where O, started to become depleted. SCOC was determined 
from Fick’s first law of diffusion. 


Microbiology sampling 

Noduleand sediment samples for microbial community analyses were 
collected from the 5D experimental chambers. Approximately 30 g 
of sediment from each of the 0-2 cm and 2-5 cm horizons and 50 g 
of intact nodules were placed in separate sterile Whirlpak bags with 
a pre-sterilized spatula and then transferred to a -80 °C freezer. DNA 
from approximately 10 g of nodules and 250 mg of sediment were 
extracted using the Qiagen PowerMax soil and PowerSoil extraction 
kits, respectively. Extracted DNA was then shipped on dry ice to Lara- 
gen Inc. and sequenced using a proprietary in-house method. The V4 
region of 16 SrRNA genes were amplified using the Earth Microbiome 
Project protocol” with the 515 F (5’-GITGYCAGCMGCCGCGGTAA”) 
and 806 R (5’-GGACTACNVGGGTWTCTAAT~) primers. Raw fastq 
files were processed using a custom pipeline (https://github.com/ 
Boston-University-Microbiome-Initiative/BU16s) built with QIIME 
2020.2 (https://www.nature.com/articles/s41587-019-0209-9). Adap- 
tor sequences were removed using cutadapt (https://doi.org/10.14806/ 
ej.17.1.200), read truncation positions were determined by mineer 
(more below), amplicon sequence variants (ASVs) were generated 
using dada2 (trunc-len-r”°) (https://doi.org/10.1038/nmeth.3869) 
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and ASVs were clustered to 99% identity with the SILVA 132 database 
(https://academic.oup.com/nar/article/42/D1/D643/1061236) using 
the vsearch cluster-features-closed-reference (https://doi.org/10.7717/ 
peerj.2584). Due to drops in sequencing quality, all reverse reads were 
truncated by 49 bases (from a length of 301 to 252) as determined by 
minERR, an algorithm for determining optimal sequence length based 
on sequence quality scores (https://github.com/michaelsilverstein/ 
mineer). Family- and genus-level abundance was computed by sum- 
ming the relative abundance of all ASVs with the same family/genus 
classification within each sample. Spearman correlations were then 
computed between family- and genus-level abundance and observed 
optode-derived total O, changes. Sequences have been archived at 
National Centre for Biotechnology Information GenBank under the 
Bioproject ID PRJNA1117483. 


Polymetallic nodule surface area measurements 

Photographs of the surface sediment and nodules in the chambers 
were imported into ImageJ. The outline of each nodule in each cham- 
ber photograph was then traced and the surface area of the nodule 
automatically calculated in Image J (assuming each surface nodule 
was flat in shape) and logged as an ImageJ file before being exported 
and saved as an Excel file. 


Radiolysis O, production estimates 

To estimate the potential radiolytic O, production, published concen- 
trations of *°U, 7*U,?”Th, *°K (refs. 22-26) in seawater were used (Sup- 
plementary Table 1). For nodules, 7°U, 7*>U and 2”Th isotopes of three 
nodules from chamber experiments from the 5D cruise were measured 
by Multicollector-Inductively Coupled Plasma Mass Spectrometer 
using previously described methods” *’ and averaged; “°K values were 
derived from the literature”. Nodule and seawater contributions were 
calculated using a kinetic model developed by ref. 9 that incorpo- 
rates 32 reactions (equation (1) in ref. 30). The nodule boundary layer 
was assumed to be fully integrated with the seawater, surpassing the 
respective ~23 to ~452 um stopping power distance of alpha and beta 
particles used to model geologic materials”. Sediment radiolytic O, 
was calculated as half of the previously quantified H, production rates 
in equatorial Pacific subsurface sediment”, given the stoichiometry of 
water’s radiolytic decomposition (an equivalency that probably offers 
an overestimate of derived O,). Contributions from these three com- 
ponents (nodules, sediment and seawater) were scaled by the benthic 
chamber’s size and contents to produce an estimate of 0.18 umol lof 
O, generated over 48 h according to the following expression. 


(02), = (Qiz x Ea X G(O2) x Mon X Aiz X 107°) x (1 — e72) 


Here (O,), is the mass (kg) of O, produced over a given time t (yr), Qz 
is the mass (g) of the isotope, E, is the average energy (eV) released 
from the decay of one atom; G(O,) is the radiation chemical yield of 
molecules per 100 eV of the radiation energy; M,, is the O, molecular 
mass (g), A is the isotope atomic mass (g) and Ais the isotope-specific 
decay constant (y*). The overall (O,), value summed the contributions 
from ”®U,’””U, Th and “°K across water, nodule and sediment sources. 


Electrochemistry measurements 

Voltage potentials were measured using a Keithley DMM6500 digital 
multimeter on nodules previously collected by coring in the UK1, 
NORI-D and BGR license areas. Nodules were initially immersed for 
seven days in Instant Ocean artificial seawater (salinity 35). Tomeasure 
the potentials, two electrodes (platinum wire, 99.9% purity) were first 
washed in perchloric acid, rinsed in Milli-Q water and dried before 
being attached to alligator clamps attached to the multimeter. The 
platinum wires were then immersed in Instant Ocean artificial seawater 
inaglass petridish to measure background voltages (0.003 + 0.001 V, 
SE, n=17) until stable. Once stable, a nodule was placed in the petri 


dish and the platinum probes placed on the nodule at random loca- 
tions, ensuring contact in one of two ways. We either carefully drilled 
ahole into some nodules so one platinum wire could be fixed inside it 
while the second platinum wire was firmly pressed against the nodule 
surface using a clamp. Alternatively, the platinum wires were pressed 
firmly against two different spots on the nodule surface and held in 
place using a clamp. Voltages were then recorded for 1-2 min until the 
signal was stable. This procedure was repeated up to 20 times in differ- 
ent randomly selected regions of the nodules depending on their size. 
Measurements were undertaken on 12 nodules at 21 °C (n =153) anda 
single control rock composed of metamorphosed carbonate (n = 10). 
Two nodules from UK1 were also retested after being cooled to 5 °C 
(n= 18) by placing them in Instant Ocean water in a refrigerator over- 
night. Voltage potentials (n = 20) between two nodules were measured 
using four nodules collected from UK1. Potentials measured during 
each measurement were averaged and corrected for the background 
seawater voltage measured using only Instant Ocean seawater in the 
absence ofa nodule. Measured resistances inside some of the nodules 
that were broken up were in the KO to 100s of KQ range, though it is 
unclear if these resistivities change at the nano- or microscale requir- 
ing further investigation. 


Geochemistry modelling 

The chemical stability and solubility of manganese (IV) oxide (birnes- 
site) to dissolved Mn” as a function of pH and O, activity was modelled 
using the Geochemist Workbench Professional (version 12) software, 
with the in-built and internally consistent THERMO database. The 
conditions used for generating the phase diagram (Extended Data 
Fig. 5) represent bottom seawater as measured in the eastern CCZ with 
a temperature of 1.6 °C and chlorine and manganese concentrations 
of 0.55 M Cland 2e™® M Mn, respectively. 


Ex situ core incubations 

Opportunistic ex situ experiments were undertaken during the 5D 
cruise using sediment cores retrieved by a multi-corer from the CTA 
area (Extended Data Fig. 1). Immediately after the multi-corer arrived 
back at the surface, cores were removed and transferred to a cold lab 
held at in situ temperature. The cores were then exposed to the fol- 
lowing five treatments (administered using a 60 ml syringe), which 
included (1) Na,HCO, (0.3 uM final concentration, n = 3), (2) NH,Cl 
(10 uM final concentration, n = 3) and (3) NH,ClI (50 uM final concen- 
tration, n = 3), (4) 0.3 uM Na,HCO, + 10 uM NH, CI (final concentration, 
n=3) and (5) HgCl, (1.1 pM final concentration, n = 3). No-injection 
controls (n = 3) were also performed and separate core experiments 
in which four nodules were incubated for 48 h by themselves with no 
additions. After addition, the water phase of each core was stirred and 
a50-ml sample of top water was taken for microWinkler analysis (as 
above). Stoppers were then placed on the top of the cores, ensuring 
no air bubbles were present. The stoppers were secured tightly and the 
cores fully submerged in a large bucket containing 0.45-pm-filtered, 
cold, surface seawater (salinity 35). The bucket was covered with five 
black plastic bags and secured in the cold room with the lights turned 
off. After 48 h, the cores were removed from the bucket, and the cores 
were inspected for the presence of air bubbles. Only one core, a HgCl, 
treatment, hada gas bubble beneath the bung, which was rejected from 
further analysis, leaving n = 2 for this treatment. The other cores were 
then re-sampled for dissolved O, and analysed as before. Core-specific 
water volume measurements were used together with the change in O, 
concentration to calculate the total net O, change per core. 

To determine if our ex situ DOP detection was affected by intrusion 
of O, from the atmosphere into the core tube, two controls were per- 
formed: a shipboard test with an O, microprofiler and alab-based test 
using the Winkler method. Shipboard, aclean core tube was filled with 
Milli-Q water and sparged with N, for 10 min before beginning the test. A 
Metrohm 8663 Multimeter was inserted through a predrilled hole inthe 
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rubber stopper, allowing for O, concentration to be recorded every 5s. 
Anincrease from 39 to 69 umol I was observed over -5 h, correspond- 
ing toarate of 0.14 mmol m°? dor 4% of the3.5 mmol md‘ mean net 
DOP measured in the ex situ experiments. Back in the home laboratory, 
three of the original core tubes were filled with 4 °C, 0.2-m-filtered 
artificial seawater (salinity 35) and sparged with N, for 8 min through 
afiltered pipette tip to achieve an initial dissolved O, concentration of 
-100 pmol I (for example, the approximate starting O, concentrations 
for the shipboard experiments). The tubes were sealed with rubber 
stoppers and electrical tape, being careful to avoid bubble forma- 
tion. They were then submerged in a 32-gallon plastic garbage can of 
unfiltered seawater (O, concentration: 228.12 umol I”) in a dark cold 
room (8 °C) for 48 h. After 48 h, the tubes were quickly unsealed and 
analysed one at a time to prevent additional O, dissolution from the 
air. A 50-ml sterile syringe was used to slowly collect 10 ml of seawater 
from the centre of the core tube, being sure to avoid bubble entrain- 
ment into the syringe. The sample was carefully expelled into a10-ml 
reaction vial and fixed using the adjusted values for a10-ml sample 
according to a volume-scaled Winkler titration protocol” and the 
reagents from the LaMotte Dissolved Oxygen Test Kit. The fixation 
of each collected sample was done in less than 2 min in a fume hood. 
Dissolved O, increased by 0.11 mmol m”? d” during the 48 h, which 
corresponds to between 3.2% of the mean net DOP rate observed in 
the ex situ experiments (3.5 mmol m”? d”). Both of our control experi- 
ments provide high confidence that the diffusion of external O, into 
the core tubes did not cause the O, production measured in the ex situ 
core incubations. 


Calculations to quantify intrusion of O, from the 
polyoxymethylene chambers and lids 

Oxygen intrusion was estimated from Stephens™ who calculated that 
20.66 pmol I of O, could diffuse out of 428 cm? of polyoxymethylene 
plastic when immersed for 48 h in hypoxic water (O, diffusion rate: 
0.02 pmol O, cm” d”). To determine the total area of plastic that would 
be available for diffusion (869-1,584 cm°), we added the surface area 
of the lid to the surface area of the four walls that would be exposed 
at the seafloor (based on the depth of the water phase—above). The 
minimum and maximum areas available for diffusion were multiplied 
by 0.02 umol O, cm d” to estimate that 41.9-76.5 umol O, I” would 
diffuse out of the polyoxymethylene chamber walls and lidin48 h under 
hypoxic conditions. Thus, we are highly confident that O, leakage from 
the plastic chambers could not replicate the high O, concentration seen 
in some of our oxygenated experiments (Fig. 1). 


Data availability 

Source data are provided with this paper. These data are also available 
via Dryad at https://doi.org/10.5061/dryad.tdzO8kq6w (ref. 35), and 
geological samples were exported in accordance with relevant permits. 
The nucleotide sequences generated by metagenome sequencing have 
been deposited in the National Centre for Biotechnology Information 
database under BioProject ID PRJNA1117483. 
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Extended Data Fig. 1| Benthic chamber lander and multi-corer deployment (b-d) of NORI-D in the central abyssal Pacific. The deployment location for 
locations across the CCZ. Benthic chamber lander (BCL) locations in APEIs the multi-corer (MUC) that sampled sediments for the ex situ experiments 
1,4, and 7 (western CCZ), UK1 and OMS and NORI-D (stars) (a) and both conducted during the 5D cruise is also shown (c). 
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dead-algal biomass during expeditions 5D (A), 5E (E), and 7A (G), DIC + NH," 
during expedition 5D (B), 0.45-um filtered seawater during expedition 5D (C), 
and control (no injection) during expeditions 5D (D), 5E (F), and 7A (H). 
Each datapoint is the mean of two Winkler measurements. 


Extended Data Fig. 2| Oxygen concentrations measured from water 
samples by Winkler titration during the NORI-D benthic chamber lander 
experiments. Mean O, concentration (umol L”) measured by micro-Winkler 
analysis conducted on water samples that were collected periodically from the 
chambers through time (hr) under different treatments. The treatments were 
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Extended Data Fig. 3 | Oxygen optode concentrations measured during 
benthic chamber lander experiments in the UK1 and OMS license areas 
and APEIs1, 4, and 7. Oxygen optode readings through time (hr) from 36-hour 
abyssal (4037-5216m) in-situ benthic chamber lander experiments conducted 
in the UK1 and OMS license areas in 2015 and APEIs 1, 4, and 7 in the western 
CCZ inJune 2018. The experiments that were conducted were identical to those 
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carried out at NORI-D. The O, concentrations recorded by the optodes in the 2015 
and 2018 experiments were derived from factory calibrations undertaken 4-6 
months prior to the expeditions as in-situ temperature could not be replicated 
onboard during the optode calibration process. As such, only relative changes in 
O concentrations can be interpreted. 
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Extended Data Fig. 4 | Bar chart showing total net O, production in ex situ difference in O, concentration of the water phase overlying the sediment 
sediment cores. Mean total net O, production (umol O, core”) measured on between t =0 hours and 48 hours accounting for the core volume. Error bars refer 
sediment cores (n=1-3) exposed to a variety of treatments during 48-hr ex situ to +1standard deviation. Individual fluxes from the ex-situ incubations are also 
incubations that were carried out on the ship at in-situ temperature and in shownas data points overlying the bars. 


the dark during the 5D cruise. Oxygen production was determined from the 
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Extended Data Fig. 5 | Phase stability and solubility of birnessite in seawater 


as a function of O, activity and pH. The phase stability and solubility of 


birnessite (manganese [IV] oxide) in seawater as a function of O, activity and pH 
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Extended Data Table 1| In-situ benthic chamber lander deployment locations and depths in NORI-D 


Cruise Date Lander deployment Area Station Depth (m) 
5D May-June 2021 AKS268 CTA STM-001 4285 
5D May-June 2021 AKS271 CTA STM-001 4284 
5D May-June 2021 AKS273 CTA STM-014 4306 
5D May-June 2021 AKS276 CTA STM-014 4306 
5D May-June 2021 AKS279 CTA STM-007 4280 
5D May-June 2021 AKS282 PRZ SPR-033 4245 
5D May-June 2021 AKS286 PRZ SPR-041 4127 
5E November-December 2021 AKS313 CTA STM-014 4304 
5E November-December 2021 AKS316 CTA STM-001 4285 
5E November-December 2021 AKS318 CTA STM-007 4277 
5E November-December 2021 AKS321 CTA STM-001 4285 
5E November-December 2021 AKS328 PRZ SPR-033 4243 
TA ‘August-September 2022 AKS329 CTA TF-021 4289 
7A August-September 2022 AKS331 CTA STM-001 4286 
7A August-September 2022 AKS334 CTA TF-028 4278 
TA ‘August-September 2022 AKS336 CTA TF-021 4271 


In-situ benthic chamber lander deployment locations and depths from the 3 cruises to the NORI-D license area. Deployment areas are shown in Extended Data Fig. 1 and are defined as the 
Collector Test Area (CTA) and Preservation Reference Zone (PRZ). 
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Extended Data Table 2 | Total net oxygen change measured by O, optodes during the in-situ benthic chamber lander 
experiments at NORI-D 


Volume pa Total Oz Weight (8) TER Net_ 
Lander Chamber | Optode Treatment of water optodes production of nodules (mmol O; production/ 
deployment Ps phase p! in 48hrs (sediment noa? | respiration 
sensor logged Oz P m? d’) a 
(L) data. (pmol 02) horizon) dominated 
AKS268 2 B Dead-algal biomass NA 0-47 = = = Production 
AKS268 3 A Dead-algal biomass 3.812 0-47 1008 782 (0-Sem) 10.5 Production 
AKS271 3 A Dead-algal biomass NA 0-47 - - - Production 
AKS273 3 A Dead-algal biomass 4.598 0-47 1545 653 (0-Sem) 16.1 Production 
AKS276 2 B DIC+NHs NA 0-47 > = = Production 
AKS276 3 A DIC+NH; 3.933 0-47 671 390 (0-Sem) 61 Production 
‘AKS279 3 A Filtered seawater 4.659 0-47 1639 410(0-Sem) 16.5 Production 
AKS282 3 A Dead-algal biomass 4.598 0-47 246 - 17 Production 
AKS286 2 B DIC+NH«" NA 0-45.48 - - - Production 
AKS286 F A DIC+NH«' 6.050 0-47 329 466(0-Scm) 24 Production 
AKS313 1 B Control (no injection) NA 0-47 - - - Production 
AKS313 2 A Control (no injection) NA 0-47 > = = Production 
AKS316 J B Control (no injection} 2.118 0-47 639 512 (0-2cm) 59 Production 
AKS316 2 A Dead-algal biomass 2.420 0-47 1526 662 (0-2cm) 15.6 Production 
AKS318 1 B Control (no injection) NA 0-47 = = = Respiration 
AKS318* 3 A Control (no injection) 2.783 0.23 -47 - 596 (0-2cm) - Production 
AKS321 F B Control (no injection) NA 0-47 a a a Respiration 
AKS321* 3 A Dead-algal biomass 2.783 133-47 - 863 (0-Sem) = Production 
AKS328 2 B Dead-algal biomass 3.691 0-47 257 630 (0-Scm) 2.0 Production 
AKS328 3 A Dead-algal biomass 3.812 0-47 683 556 (0-Sem) 64 Production 
AKS329 1 A Control (no injection) 3.872 0-47 904 755 (0-5cm) 95 Production 
AKS329 Z B Control (no injection) 3.872 0-47 774 704 (0-Sem) 8.1 Production 
2 A Dead-algal biomass 2.723 0.06 -46.45 > 704 (0-Sem) = Production 
3 B Dead-algal biomass 3.207 007-47 = 608 (0-Sem) = Production 
2 A Dead-algal biomass 4.175 0-47 1710 667 (0-Sem) 18.0 Production 
AKS334 3 B Control (no injection) 5.143 0-47 652 621 (0-Sem) 5.7 Production 
AKS336 3 B Control (no injection) 4.477 0-47 1408 622 (0-5cm) 14.4 Production 
Mean + SE 933+ 132 93415 
Total net oxygen change measured by O, optodes during the in-situ benthic chamber lander experiments at NORI-D. Stars (*) denote the incubations where the optodes did not start to log 
O, concentrations immediately or stopped logging before 47 hr meaning the total O, change could not be determined (designated by a dash [-]). The total net O, change could also not be 


determined in a sub-set of experiments (designated by a dash [-]) because the chambers failed to seal sediments at the end of the experiment, which meant that the volume of water in the 
chamber could not be determined (designated NA) once the lander arrived back on deck. DOP flux rates were determined from the change in O, concentration from the point when the 
optodes started logging to t = 47-hr. The determination of net production/ respiration dominated experiments was based on the O, optode profiles seen in Fig. 1. 
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Extended Data Table 3 | Theoretical diffusion times for thin versus thick-walled bubbles at the seafloor 


Time (sec) required for Time (sec) required for 
diffusion assuming a diffusion assuming a 
thin-walled (10 nm) thick-walled (10000 

bubb 

AKS268-Ch3 0.012 1.226 

AKS273-Ch3 1.411 

AKS276-Ch3 1.068 

AKS279-Ch3 1.442 

AKS282-Ch3 0.768 


Incubation 


le 

AKS286-Ch3 0.854 
AKS316-Ch1 1.053 
AKS316-Ch3 1.407 
AKS328-Ch?2 

AKS328-Ch3 0.01 1. 
AKS329-Chl 0.012 1.182 
AKS329-Ch2 0.01 1.122 
AKS334-Ch2 0.015 1.463 
AKS334-Ch3 0.01 1.061 
AKS336-Ch3 0.014 1.372 


0.012 + 0.001 1.153 + 0.061 


Theoretical diffusion times (seconds) for thin versus thick-walled bubbles at the seafloor. The theoretical time for O, to diffuse from a trapped air bubble at the seafloor was estimated 

by calculating the size of air bubble that would be required at the surface to produce the O, increase seen in each chamber. This was estimated from the difference between the max O, 
concentration value recorded by the optode and the initial O, optode reading in the chamber, the volume (L) of the water phase in the chamber and assuming air was comprised of 21% O,. 
The size of the bubble at the seafloor was then computed using Boyle’s law and the in-situ pressure calculated from the deployment depth. The O, concentration gradient (dC) between the 
bubble and the water phase of the chamber was determined from the concentration of O, in the theoretical bubble and the initial concentration of O, in the chamber. The diffusion distance 
(dZ) was set at 10 and 10,000 nm to calculate the diffusion times across a thin versus thick-walled bubble. The diffusion coefficient for O, (1.064 x 10° cm? s”) was calculated from the O, 
diffusion coefficient based on an in-situ salinity of 35 and temperature of 1.6°C. Fick’s first law of diffusion was then used to calculate the diffusion time based on the diffusion coefficient, 


dZ and dC. 


i=} 


0.780 
080 
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Extended Data Table 4 | Minimum and maximum voltage potentials (mV) measured on the surface of polymetallic nodules 


Min. Nod.1 | Nod.2 | Nod.3 | Nod.4 | Nod.5 | Nod.6 | Nod.6 | Nod.7 | Nod.7 | Nod.8 | Nod.9 | Nod. Nod. Nod. Cont. 
voltage (cold) (cold) 10 11 12 
(mV) 
Record 1 0.00 0.00 0.00 0.00 0.01 0.00 58.55 77.83 0.00 0.01 0.00 0.00 78.81 0.00 0.01 
Record 2 8.97 18.45 15.48 3.03 0.00 9730 | 40.90 | 79.56 | 21.65 | 25.34 0.00 0.02 53.68 0.01 4.48 
Record 3 12.00 0.00 26.48 78 3.14 100.26 | 57.69 | 81.87 22.13 5.57 0.02 031 0.00 22.30 0.97 
Record 4 8.67 1.01 29.21 21.06 | 24.87 88.69 | 54.14 | 44.84 37.31 36.89 15.42 23.26 0.00 0.00 
Record 5 1.17 0.64 18.68 21.91 14.89 | 5268 | 60.72 | 2236 30.80 50.57 0.00 9.28 1.83 3.34 
Record 6 0.00 6.28 8.45 13.81 17.71 49.47 85.98 55.17 71.57 38.20 0.02 0.00 
Record 7 0.00 2.95 3.98 14.25 0.90 6944 | 75.57 | 75.75 72.06 | 44.82 1.10 
Record 8 14.17 0.00 16.97 14.62 15.03 66.93 0.07 78.40 | 57.64 0.01 0.00 
Record 9 5.57 0.40 0.00 0.00 9.35 33.69 81.75 59.83 8.03 1.66 
Record 10 5.00 3.52 3.29 18.15 24.79 0.04 20.69 | 65.61 27.96 3.53 
Record 11 6.14 3.39 32.39 13.09 25.84 36.46 
Record 12 1.52 5.21 25.42 16.84 25.45 36.46 
Record 13 0.48 0.16 27.85 | 26.23 0.01 61.34 
Record 14 7.12 4.09 28.61 24.92 12.71 60.85 
Record 15 1.51 0.00 16.31 0.00 51.72 
Record 16 1.59 3.09 15.70 0.00 76.51 
Record 17 0.00 0.00 22.12 4.70 78.77 
Record 18 7.83 0.00 13.64 17.90 65.17 
Record 19 6.33 0.75 20.67 24.68 49.45 
Record 20 10.06 12.24 29.94 20.49 49.45 
Max. 
voltage 
(mV) 
Record 1 10.93 52.73 | 266.80 | 52.99 | 45.69 | 102.87 | 70.85 80.11 | 116.23 | 39.27 | 219.67 | 19.36 | 98.41 | 128.40 | 52.91 
Record 2 9.81 21.46 | 22.95 4.89 0.47 99.85 42.24 | 8260 | 30.74 | 28.94 | 72200 | 112.63 | 70.40 | 407.56 72 
Record 3 13.84 12.00 | 28.36 3.01 8.44 103.48 | 59.12 | 82.80 | 27.61 6.61 485.52 | 178.04 | 952.61 | 59.50 11.15 
Record 4 12.28 3.51 31.39 | 24.88 | 28.98 | 96.82 | 58.33 57.64 38.58 39.47 | 350.78 | 65.84 600.43 2.96 
Record 5 8.99 4.87 75.42 23.15 21.01 77116 | 6541 3191 34.22 53.26 | 361.71 | 25.45 118.71 7.56 
Record 6 3.17 7.27 12.12 14.86 | 25.77 | 74.55 | 266.76 | 80.50 | 91.66 | 42.18 128.61 8.71 
Record 7 4.29 5.08 9.94 14.86 111 7773 | 102.45 | 78.95 73.12 77.14 2.85 
Record 8 25.22 2.63 20.15 18.49 16.85 69.20 | TAI 84.31 5791 11.12 0.88 
Record 9 7.45 134 9.36 8.69 10.91 37.81 82.78 65.62 11.87 2.45 
Record 10 8.38 14.61 3.81 28.73 26.73 | 113.62 28.34 | 65.91 33.35 8.70 
Record 11 173 9.19 35.00 14.85 28.08 39.46 
Record 12 3.17 7.02 32.99 18.21 33.05 39.46 
Record 13 2.94 5.52 32.21 31.65 12.76 63.50 
Record 14 9.83 6.43 31.86 | 26.59 21.71 61.74 
Record 15 2.45 7.04 22.45 3.49 52.38 
Record 16 6.76 6.45 21.56 11.27 79.14 
Record 17 5.54 437 24.97 9.10 81.13 
Record 18 10.19 0.95 20.30 20.38 68.18 
Record 19 7.91 2.13 24.16 27.89 52.58 
Record 20 14.32 13.57 34.40 27.74 52.58 
Minimum and maximum voltage potentials (mV) measured on the surface of the polymetallic nodules (nod.) at 21°C (nodules 1 - 12) and 5°C (nodules 6 and 7 cold), and a control (cont.) rock. 


The voltages have not been corrected for the background voltages measured using only instant water. 
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